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The exact expressions for density matrix and Wigner 
functions or quantum systems are known only in special cases 
and .practically, all or them and their references are 
described in [1-31. Correspondending Hamiltonians are 
quadratic forms of Euclidean coordinates and momenta. In 
this paper we consider the problem or one-dimensional free 
particle movement in the bounded region 0 < x < a (including 
the case a= ®). For this problem the solutions of Sc hr odinger 
equation are well known: 

v' n =(2/a) 1 '' 2 sin(nnx''a), E n =(rrnh) 2 -'(2ma 2 >. n-1.2.3.... (1) 

Then the equilibrium density matrix can be calculated by 
formula 

p(x.x’ v'hWv'^x* ) exp(-/?E n ) (2) 

Indroducing the expression (1) to (2) and mading some 
simple transformations we obtain two series, each of which 
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is, practical liy the definition or teta-function [41 



Hie replacement transforms the expression (4) into 

the propoaator or Schrodinger equation for the particle in a 
box obtained earlier by various methods in [3,5,71. 

Evidently in the limit ( ^a ej - <» (i.e. low temperature 

and small size of a box) the density matrix can be good 

approximated only by the first order term of the expansion 

series (2). The question is in obtaining from the exact (but 

not very obvious) formula (4) the asymptotics or the density 

7 

matrix in quasi - classical limit ft' a (high temperature and 
wide box). 

The qualitive behaviour of the probability density 
p(x,x,ft) in this case is clear from physical consideration. 
It must be almoust constant at all points inside the box 
except very small region near the wall corresponding to de 
Brogie wave length. In this region the density matrix must 
leads to zero- However, it is interesting to obtain this 
result from the formula (4). More over we would like to know 
the character of the deflexion uniform distribution inside 
the box caused by quantum corrections. This problem can be 
solved using the equality 141 for teta-function 


^ Vj exp (-niz^Vr )© 3 (-z/t | -l'-T ) 


Due to the fact that in our case the parameter r is pure 
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complex and restricting by the first term of expansion 
series (3) of the function |-i^r) when r-*o we obtain 

the following formula describing the quasi -cl assi cal 
behaviour of the density matrix 


p(X. X’ . ft)- (m''2TT fib! 


r 

n * O 
7 


exp [ -m(x-x’ (1+ 


2exp [ -ma^h^lch [-2m(x-x’ )/&**] )-exp [ -m(x+x’ )*v2/?h*3 (1 


2exp [ -ma 2 -Y5^ 2 lch [-2m(x+x’ )| 


(6) 


This formula is correct in the region |x + x’ a (i-e. at 
the left half of the box ). For the points outside of this 
region one have to use the properties folloing from (3) 
and (4) 

p(x,x* )=p(x’ ,x). p(a-x,a-x’ )=p(x,x’ ) 

For the diagonal elements of probability density we have 
following expression 


P (x, X’ = (m- 2n/?h } 


n a ,n 

£r \ ±' 


1-exp [ -2mx 2 'Y3h 2 j + 2exp i-ma^WJh^l [1- 


- exp [-Smx^N 3ch Umx^h 1) 


x<8'2. 


>> 1 ( 8 ) 


The first two terms in figure brackets describe 
the probability density of particle position in the infinite 
half space right from the wall placed to the point x=0. The 
other terms give corrections caused by the-presence of the 
second wall. Note, that this corrections don’t oscillate as 
it can be seemed from formulas (3) and (4). 

In the centre of the box the density matrix is equal to 
p (a-Z, a^2./3)=const <1-2 exp(-ma*v2^/ J ) 
and the half space case on the same distance from the 
coordinate centre we have an analogous expression but 
without two in front or exponent- The exact expression of 
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statistical sum have the form 


z (*)= — 


Li ,i h .1 

U it -1 

l 2ma J J 


/ *S \ 

i?,' 


It's quasi -cl assical expansion is 


^ * r\ 

£ ( r ?} = 8 [U\ y ' fjJi }*h ~ j (_t‘ , *io6Xp l y j fVH^^ ^2 (±0) 




From (4) one can obtained the Wigner function 


W(p,q,fS)= lo(q+jf/2.q-Jf/2,/?) exp(-ip^f/h) dr 


(11) 


Taking into consideration that the integration region is 
bounded by the interval -2q < ? < 2q under the condition 
0 < q < a/2 we have [71 


W(p,q,/?)= 



(2py/h) ©o 


- (h/ap)sin (2pq/h)© 0 

KJ 


f y 

irrh ft 2 ' 

l"a 

2ma 2 . 

r q 

inh o 2 

— 

rr 

. a 

2ma'-' „ 


( 12 ) 


but when a/2 < q < a one have to use the equality 

W (p, q, /?) = W(p,a-q,f5) 

The Wigner for a free particle in half space was exactly 
expressed by the error -function for the first time in [81 
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